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APPENDIX C 
PRELIMINARY CHANNEL STABILITY ASSESSMENT 



C.l CHANNEL STABILITY APPROACH AND METHODS 

Several approaches exist for developing design criteria and evaluating channel stability for river 
reconstruction efforts. Approaches are generally categorized as analog, empirical or analytical 
methods. Analog methods rely on field measurements from stable reference reaches to develop 
dimensionless ratios for pattern, cross section and profile characteristics that are applied to 
design reaches. Empirical methods use regime equations and regional hydraulic geometry 
relationships to predict hydraulic properties. Analytical methods employ numerical models to 
simulate hydraulic conditions for a range input variables. Due to theoretical assumptions and 
limitations that are unique to each method, an unacceptable range of answers typically results. 
Therefore, the standard channel stability assessment involves using a combination of several 
methods and relying on experience, practicality and judgment to interpret the results. 

In the absence of a standardized approach to stream channel design, the design team used a 
combination of elements from several techniques that represent the best available methods for 
developing design criteria for restoration designs. Interpreting results, measuring channel 
stability and establishing acceptable design thresholds relied on the experience and judgment of 
the designers, the application of standard hydraulic principles and not necessarily an established 
or accepted set of design criteria. As described in Section 3, analog, empirical, and analytical 
methods provided the basis for developing a range of design channel dimensions, and were used 
to predict the most probable cross section, plan form and profile dimensions for the CFR and 
BFR. 

Final channel dimensions were established through a trial and error process that included several 
iterations of channel stability analyses. The methods used to complete the channel stability 
analysis focused on refining the cross section, plan form and profile dimensions until acceptable 
values were observed for each hydraulic parameter considered to influence the dynamic 
equilibrium of river hydraulics. Acceptable values were assumed to be a range of values that 
best achieved dynamic equilibrium among hydraulic parameters, anticipated sediment size and 
anticipated sediment supply. Emphasis was placed on the results of analog methods for the 
selection of cross section, plan form and profile dimensions with considerable consideration 
given to sediment supply and sediment size. A range of +/-15% around a mean value for a 
specific hydraulic parameter was considered to be the guideline for the variation that a newly 
constructed river can accommodate before dynamic equilibrium is disrupted and instability is 
triggered. Although a statistical analysis was not employed, this guideline, along with additional 
consideration of reference reach conditions, was employed as threshold design criteria for 
channel stability and maintaining dynamic equilibrium. For reporting purposes, only the final 
iteration of analyses was reported in Appendix C - Preliminary Channel Stability Assessment. 

In terms of the analyses in Appendix C, channel stability represents a condition where several 
hydraulic variables are balanced to achieve a state of dynamic equilibrium that approximates 
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reference conditions, satisfies traditional hydraulic design principles and considers results from 
the best available regime equations. Channel stability implies that although the channel pattern 
may change over time, the channel's cross section area and slope remain consistent. Under 
stable conditions, rates of erosion and deposition are approximately balanced as the channel 
pattern changes. 

The scale of this project along with the public health risks and implications of failure on 
surrounding infrastructure resulted in an intensive data collection effort that allowed the use of 
several methods to evaluate channel stability and predict stable, post-dam-removal river and 
floodplain conditions. As discussed in Appendix B, significant effort was devoted to data 
collection on selected reference reaches on the CFR and BFR both upstream and downstream of 
the restoration project area. These reaches are assumed to be stable and functioning in a state of 
dynamic equilibrium. In addition to displaying stability, the reference reaches possess the 
desired biological environment for aquatic wildlife and riparian vegetation species. For these 
reasons, emphasis was placed on methods that derived results from reference reach data. Figure 
C-1 identifies the location of the selected reference reaches in relation to the restoration project 
area. Sheet 1-7 in Appendix I identifies the locations of reference reaches in CFR3. 

The BFR Gage, CFR3 Middle (CFR3-B, see Figure 1-1) and CFR Bandmann Flats reference 
sites were identified as best examples of the most probable channel conditions for the BFRl 
reach, CFR2/CFR3 reach and CFRl reach, respectively. To validate the assumption that the 
selected reference reaches are stable, a comparison of hydraulic properties predicted by methods 
using analog, empirical and analytical approaches was completed. Results were used to establish 
a range of design criteria for proposed cross section and profile dimensions for the BFR and 
CFR. 

The following assessments focus on validating the feasibility of the conceptual plan. As such, 
assessments were based on conceptual design layouts. Specific hydraulic properties evaluated 
included slope, velocity, shear stress, width-depth ratio, and sediment transport (power). 
Hydraulic properties were computed for a range of discharges from bankfull discharge to 100- 
year discharge. Additional consideration was given to effects of ice floes and ice jams on 
channel stability. An assessment of the selected channel dimensions over a range of discharges 
was completed using HECRAS and is presented at the end of this appendix. Hydraulic modeling 
efforts undertaken for proposed conditions used simplified channel cross sections and floodplain 
geometry to estimate potential design conditions. Hydraulic parameters derived from the analysis 
were used as input for calculations presented in the channel stability assessment. Further 
refinement of these assessments is expected during the final design phase. 
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C.2 STABLE SLOPE ANALYSIS 

Formation of equilibrium channel slope is influenced by valley longitudinal slope, extent of 
lateral channel confinement within the valley and substrate composition. Generally, methods for 
estimating stable channel slope predict slopes that correspond to a threshold hydraulic condition 
for the initiation of movement for a specific bed material size class. As such, available equations 
rely on input data for bed composition, mean depth and dominant discharge. 

C.2.1 Methods 

A stable slope analysis was completed for selected reference reaches. Estimated slopes were 
compared with actual average slopes and used to validate proposed channel slopes. For the 
purpose of this analysis, it was assumed that available bed material in the restoration project area 
will resemble that found in the corresponding reference reaches. The following five methods 
were selected to evaluate critical slope. 



Schoklitsch(1932) 
S = K(D5o/Qf'' 
S = Slope (ft/ft) 
K = Constant, 0.00174 



D50 = Particle size for which 50% is finer (mm) 
Q = Dominant discharge (cfs) 



Meyer-Peter and Muller (1948) 

S = [K (Q/Qb)(n,/(D9o''')'''D5o]/d 

S = Slope (ft/ft) 

K = Constant, 0.19 

ns = Manning's roughness for bed 



Q = Total discharge (cfs) 

Qb = Discharge over bed of river (cfs) 

D90 = Particle size for which 90% is finer (mm) 



D50 = Particle size for which 50% is finer (mm) 



MacB room (1981) 
S = (KWbfDso/Q)""' 
S = Slope (ft/ft) 
K = Constant, 0.00021 
Wbf = Bankfull width (ft) 



D50 = Particle size for which 50% is fmer(mm) 
Q = Dominant discharge (cfs) 



Rosgen(2001) 

S = T Ys Di/DS = Slope (ft/ft) 

c 

D = Mean depth (ft) 

Ys = Specific weight of sediment (lbs/ft^) 



T = critical dimensionless shear stress 

c 

Di = Pavement Dioo or bar sample Dioo (ft) 



Thorne and Hey (1997) 

S = 0.087Q-°'4^D5o-"°'D84°-^V-'" 
S = Slope (m/m) 
Q = Bankfull discharge (cms) 
Qs = Sediment discharge (kg/m) 



D50 = Particle size for which 50% is finer (mm) 
D84 = Particle size for which 84% is finer (mm) 
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US Army Corps of Engineers (1994) Nomograph 
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Figure C-2. Guide to slope-discharge relationships for erodible channels. 

As with other empirical methods, these equations have limitations of applicability. Many have 
been developed for the purpose of designing armored flood control channels. Also, some 
equations incorporate factors of safety to ensure that bed movement does not occur. In addition, 
some equations have been calibrated to uniform particle sizes and sand bed channels. For these 
reasons, a wide range of results can be expected. 

C.2.2 Results 

Table C-1 presents the results of the analysis. 
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Table C-1. Comparison of estimated stable slope and average slope (ft/ft) for selected reference reaches 
and proposed reaches. 


Critical Slope 
Method 


CFR3 
Middle C4 


CFRl 
Bandmann F3/1 


BFRl 
Bonner B3 


BFR 

Ovando 

C4 


BFR 
Ovando B3 


BFR 
Ovando F4 


Schoklitsch 


0.0033 


0.0028 


0.0049 


0.0029 


0.0037 


0.0027 


Meyer-Peter 
and Muller 


0.0047 


0.0044 


0.0098 


0.0069 


0.0092 


0.0067 


US Army Corps 
of Engineers 


0.0027 


0.0020 


0.0055 


0.0030 


0.0035 


0.00250 


MacBroom 


0.0033 


0.0028 


0.0049 


0.0029 


0.0037 


0.0027 


Rosgen 


0.0035 


0.0033 


0.0063 


0.0048 


0.0056 


0.0041 


Thorn and Hey 


0.0043 


0.0044 


0.0064 


0.0043 


0.0045 


0.0041 


Average of 
Methods 


0.0040 


0.0034 


0.0061 


0.0039 


0.0047 


0.0042 


Actual Average 
Slope 


0.0028 


0.0019 


0.0032 


0.0024 


0.0025 


0.0024 


Selected Design 
Slope 


0.0031- 
0.0043 


0.0030 - 0.0032 


0.0025 - 
0.0040 


N/A 


N/A 


N/A 



C.2.3 Conclusions 

Results indicate that actual average slopes are less than estimated stable slopes for nearly all 
methods and all reaches. Since the selected reference reaches are stable and do not appear to be 
aggrading, it can be concluded that a range of acceptable slopes is available for the proposed 
design reaches. 

Existing project constraints and required tie-in points throughout the restoration project area 
dictated the values for selected design slopes in Table C-1. Despite these limitations, most 
selected design slopes lie within the acceptable range for stable slope. For those values less than 
the acceptable range, it may be necessary to adjust other hydraulic variables such as width-depth 
ratio so that the proposed channel does not experience aggradation. For those values greater than 
the acceptable range, it may be necessary to increase the frequency of grade control structures or 
increase the width-depth ratio. 

C.3 VELOCITY ASSESSMENT 

Channel velocity is a function of channel slope, channel geometry and channel roughness. 
Available methods for estimating critical velocity predict mean velocities that correspond to a 
threshold hydraulic condition for the initiation of movement for a specific bed material size 
class. As such, available equations rely on input data for bed composition, mean depth and 
slope. 
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C.3.1 Methods 

A critical velocity analysis was completed for selected reference reaches. A combination of 
analog, empirical and analytical methods was used. Hydraulic models used included WinXSPro 
(WEST Consultants, Inc. 1998) and HEC-RAS (US Army Corps of Engineers 2004). Both 
models were calibrated by using field-surveyed channel geometry, bankfull indicators and 
bankfull discharges to predict Manning's roughness values. Table C-2 presents a summary of 
the Manning's roughness values used for the analysis. 



Table C-2. Summary of estimated Manning's roughness values for selected reference reaches and 
proposed reaches at bankfull discharge. 


Method 


CFR3 
Middle C4 


CFRl 
Bandmann 

F3/1 


BFRl 
Bonner B3 


BFR 

Ovando 

C4 


BFR 
Ovando B3 


BFR 
Ovando F4 


HEC-RAS 


0.031 
(0.025- 
0.036) 


0.032 
(0.027- 
0.037) 


0.039 
(0.038- 
0.044) 


0.032 
(0.028- 
0.036) 


0.036 
(0.034- 
0.038) 


0.035 
(0.032- 
0.038) 


Hey 


0.035 
(0.034- 
0.037) 


0.037 
(0.036- 
0.038) 


0.043 
(0.042- 
0.044) 


0.038 
(0.037- 
0.039) 


0.046 
(0.042- 
0.049) 


0.037 
(0.037- 
0.038) 


WinXSPro 


0.029 
(0.025- 
0.036) 


0.032 
(0.027- 
0.037) 


0.035 
(0.032- 
0.039) 


0.032 
(0.028- 
0.036) 


0.036 
(0.034- 
0.038) 


0.035 
(0.032- 
0.038) 


Average of 
Methods 


0.032 


0.034 


0.039 


0.036 


0.039 


0.036 


Selected Design 
Roughness 


0.032-0.035 


0.040 


0.037 


N/A 


N/A 


N/A 



In addition to estimated roughness values presented in Table C-2, bankfull roughness values 
were calculated for the BFR Bonner gage and CFR Missoula gage. WinXSPro was used to 
process field surveyed bankfull indicators and hydraulic geometry, and calibrate roughness 
values. Calculated values were 0.039 for the BFR Bonner gage and 0.030 for the CFR Missoula 
gage. 

The following five methods were selected to evaluate velocity. 

Isbash (1936) 

Dso = Particle size for which 50% is finer (ft) 



V=10.7(D5o)^'^d^''^ 



V = Permissible velocity (ft/s) 
d = Mean depth (ft) 

Mavis and Laushey (1948) 
De=1.88V^' 

V = Mean velocity (ft/s) Dc = Armor size (mm) 
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Hey (1979) 

Modified Darcy-Weisbach 

U = (8/f)^'^= 5.75 log(aR/3.5*D84) 

a=ll.l(R/dm)"*'-^^^ 

u* = (9.8RS)^°^ 

f = Darcy-Weisbach resistance coefficient 

a= Function of channel cross-sectional shape varying between 11.1 and 13.46 

R = Hydraulic radius (m) 

D84= Particle size for which 84% is finer (m) 

dm = Maximum depth at the section (m) 

U = Mean velocity (m/s) 

u* = Mean shear velocity (m/s) 



Manning (1891) 
V=[1.49R^'^S^'V 
V = Mean velocity (ft/s) 
S = Slope (ft/ft) 



R = Hydraulic radius (ft) 

n = Manning's resistance coefficient 



Estimated critical velocities were used to validate proposed channel velocities and compute 
proposed channel cross sectional area. For the purpose of this analysis, it was assumed that 
available bed material in the restoration project area will resemble that found in the 
corresponding reference reaches. As with other empirical methods, these equations have 
limitations of applicability. Many have been developed for the purpose of designing armored 
flood control channels. Also, some equations incorporate factors of safety to ensure that bed 
movement does not occur. In addition, some equations have been calibrated to uniform particle 
sizes and sand bed channels. For these reasons, a wide range of results can be expected. 

C.3.2 Results 

Table C-3 presents the results of the analysis. 

Table C-3. Comparison of critical velocity and calculated velocity (ft/s) for selected reference reaches 
and proposed reaches. 



Velocity 


CFR3 


CFRl 


BFRl 


BFR 

Ovando 

C4 


BFR 


BFR 


Method 


Middle C4 


Bandmann F3/1 


Bonner B3 


Ovando B3 


Ovando F4 


Isbash 


5.1-7.1 


5.8-9.0 


8.0- 10.6 


5.6-8.0 


6.3-9.6 


5.2-7.9 


Mavis and 
Laushey 


5.2-7.3 


5.9-9.2 


8.3- 10.9 


5.7- 8.2 


6.4-9.9 


5.4-8.1 


Hey 


5.1 


6.0 


5.4 


4.0 


3.6 


4.7 


Manning (Qbf) 


5.7 


7.2 


6.6 


4.9 


4.6 


5.0 


Manning (Qbf) 


5.5 


7.1 


6.1 


5.4 


5.2 


5.3 


Manning (Q2) 


6.2 


8.1 


6.9 


6.3 


5.9 


6.2 
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Table C-3 (Continued). Comparison of critical velocity and calculated velocity (ft/s) for selected 
reference reaches and proposed reaches. 



Velocity 


CFR3 


CFRl 


BFRl 


BFR 

Ovando 

C4 


BFR 


BFR 


Method 


Middle C4 


Bandmann F3/1 


Bonner B3 


Ovando B3 


Ovando F4 


Manning (Qio) 


7.8 


9.8 


8.2 


6.8 


6.4 


7.1 


Manning (Q25) 


8.5 


10.4 


8.7 


7.2 


6.8 


7.7 


Manning (Q50) 


8.9 


10.8 


9.0 


7.8 


7.1 


8.5 


Manning (Qioo) 


9.3 


11.1 


9.2 


8.4 


7.3 


9.1 


Manning (Q500) 


10.2 


11.5 


9.7 


10.3 


7.4 


11.0 


Average of 
Methods (Qbf) 


5.3 


6.4 


6.9 


5.1 


5.2 


5.1 


Selected Design 
Velocity (Qbf) 


5.8 


7.0 


6.4 


N/A 


N/A 


N/A 



* Critical velocity range calculated for D50 and D84 particle sizes. 

In addition to velocity values presented in Table C-3, velocity values were calculated for the 
BFR Bonner gage and CFR Missoula gage. WinXSPro was used to process field surveyed 
bankfull indicators and hydraulic geometry, and calibrate roughness values. Calculated values 
were 6.6 ft/s for the BFR Bonner gage and 5.9 ft/s for the CFR Missoula gage. 

C.3.3 Conclusions 

Results yielded an acceptable range of velocities for bankfull discharge conditions. As discussed 
in C.2, existing project constraints and required tie-in points throughout the restoration project 
area dictated the values for selected design slopes. Variations in slope between reference reaches 
and project reaches account for the difference between average calculated mean velocity and 
selected design mean velocity. 

Despite these limitations, selected design velocities lie within the acceptable range for mean 
velocity. For those values less than the acceptable range, it may be necessary to adjust other 
hydraulic variables such as width-depth ratio so that the proposed channel does not experience 
aggradation. For those values greater than the acceptable range, it may be necessary to increase 
channel roughness with larger substrate or increase the width-depth ratio. 

C.4 HYDRAULIC GEOMETRY ASSESSMENT 

The fundamental concept behind assessing stable channel geometry is to establish channel 
dimensions that function in dynamic equilibrium. Cross section geometry, in particular, is 
influenced by many variables including valley morphology, channel forming discharge, slope, 
sinuosity, sediment supply and sediment composition. There have been extensive efforts directed 
toward the development of hydraulic geometry regime equations for gravel bed rivers. Most 
regime equations predict geometry as a function of channel forming discharge. However, regime 
equations typically produce a wide range of results due to the difference in experimental 
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conditions under which they were developed. As an alternative, data from stable reference 
reaches is another useful means for predicting stable channel geometry. Other methods have 
been developed to predict stable channel dimensions for flood control channels with engineered 
hnings. 

C.4.1 Methods 

A hydraulic geometry assessment was completed for the selected reference reaches and proposed 
project reaches. The assessment compared results predicted by applicable regime equations with 
measured values from selected reference reaches. Average results were used as guidance for 
developing width-depth ratios for proposed channel cross sections. Efforts focused on 
developing stable riffle cross section dimensions. The following regime equations were selected 
for the analysis: 



Williams (1986) 
W = (Qbf/0.4)''-'' 
W = Bankfull width (ft) 



0.69 



D = 0.12W' 

Qbf= Bankfull discharge (cfs) 



D = Mean depth (ft) 



-1.10 



0.125Q*'^-^V'^^- 



,0.64 



Millar (2004) 

W = 16.5Q**'™S*'"V""" D = u.izsg-- -"S-- -|i 

W = Bankfull width (m) D = Mean depth (m) 

S = Slope (m/m) D50 = Particle size for which 50% is finer (m) 

Q* = Dimensionless discharge = Q/[D^5oV(gD5o(s-l))] 

g = gravitational acceleration (m/s^) 

s = specific weight of sediment (N/m'') 

|i'= ratio of critical bank shear stress to critical bed shear stress (1.0 for no bank vegetation, 1.4 

for established bank vegetation) 



Hey and Thome (1979) 
W = 4.33Qbf° '^ 
W = Bankfull width (m) 
D = Mean depth (m) 



D 



0.22 Qbf°-^^D5o-«-^^ 



Qbf= Bankfull discharge (cms) 

D50 = Particle size for which 50% is finer (m) 



Rosgen(2001) 
D = T YsDi/S 

c 

D = Mean depth (ft) 

S = Slope (ft/ft) 

Ys = Specific weight of sediment (lbs/ft'') 



T = critical dimensionless shear stress 

c 

Di = Pavement Dioo or bar sample Dioo (ft) 



Bray (1982) 
W = 2.38Q2 



0.53 



D = 0.226 Q2' 



0.33 



W = Bankfull width (ft) Qbf = Bankfull discharge (cfs) 



D = Mean depth (ft) 
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Simons and Albertson (1963) 

W = 2.5Qbf°'' 

W = Bankfull width (m) 



tO.36 



D = 0.43W" 

Qbf= Bankfull discharge (cms) 



D = Mean depth (m) 



Lacev (1948) 
W 



2.67Qbf 

Bankfull width (ft) 



D = Qbf/[13.5V(25.4D5o)]''' 
Qbf= Bankfull discharge (cfs) 



D = Mean depth (ft) 



C.4.2 Results 



Tabled C-4 and C-5 present the results of the analysis. 



Table C-4. Comparison of width (ft) and mean depth (ft) for selected reference reaches and proposed 
reaches. 



Geometry Method 


CFR3 Middle C4 


CFRl Bandmann F3/1 


BFRl Bonner B3 




W 


D 


W/D 


W 


D 


W/D 


W 


D 


W/D 


Williams 

Millar 

Hey and Thome 


140 
154 
135 


3.6 

3.7 
5.3 


39 
41 

25 


268 
260 
244 


5.7 
5.9 
8.0 


47 
44 
31 


202 
110 
188 


4.7 
8.0 
6.1 


43 
14 
31 


Rosgen 




4.0 






7.8 






8.2 




Bray 


208 


3.7 


57 


383 


5.3 


72 


291 


4.5 


65 


Simons and Albertson 

Lacey 

Average of Methods 


82 
151 
142 


7.1 
4.9 
4.3 


11 
31 
34 


149 

272 
263 


10.9 

6.9 

6.4 


14 
40 
43 


114 
210 
186 


9.1 

5.2 
6.0 


13 
40 
36 


Average of Reference Reaches 


141 


3.1 


45 


238 


6.2 


38 


200 


5.1 


39 


Selected Average Design Value 


150 


3.7 


40 


245 


6.1 


40 


200 


4.9 


40 



Table C-5. Comparison of width (ft) and mean depth (ft) for selected reference reaches. 



Reach 


BFR Ovando C4 


BFR Ovando B3 


BFR Ovando F4 


Geometry Method 


W 


D 


W/D 


W 


D 


W/D 


W 


D 


W/D 


Wihiams 

Millar 

Hey and Thome 


147 

90 

141 


3.8 
6.2 
5.4 


39 
15 
26 


147 

78 

141 


3.8 
6.9 

5.2 


39 
11 

27 


147 
100 
141 


3.8 

5.7 

5.5 


39 
17 
26 


Rosgen 




5.8 






6.9 






5.4 




Bray 


180 


3.3 


54 


180 


3.3 


54 


180 


3.3 


54 


Simons and Albertson 

Lacey 

Average of Methods 


107 
158 
137 


4.9 
4.8 
4.9 


22 
33 
28 


107 
158 
132 


4.9 
4.7 
5.1 


22 
34 
26 


107 
158 
139 


4.9 
5.0 
4.8 


22 
32 
29 


Average of Reference Reaches 


170 


2.9 


59 


191 


3.0 


64 


183 


3.0 


61 
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C.4.3 Conclusions 

Results produced a wide range of values. As mentioned previously, differences in experimental 
conditions under which the equations were developed is a likely cause for the variability. 
However, the average of all methods when compared to actual reference reach values produced 
an acceptable range of values. When selecting design values, consideration was given to regime 
equations however, greater emphasis was placed on observed values from reference reaches. 

C.5 SHEAR STRESS ANALYSIS 

The incipient motion for a sediment particle by flowing water is a key concept in estimating 
shear stress and sediment transport. The movement of a single particle is influenced by a range 
of factors including hydraulic geometry, slope, bed composition and bed armoring. The force at 
which point a particular size particle is set in motion is termed the critical shear stress of that 
particle size. The Shields equation is used to calculate critical shear stress (Shields 1936). 

* 

T = T (y -y)D 

CCS 

* 

where x^ is the dimensionless Shields parameter for entrainment of a sediment particle of size D, 
and Yj, and y are the unit weights of sediment and water, respectively, expressed in pounds per 
cubic foot. Generally, the parameter D is taken to be Djg, the median grain size of the bed 
sediment, and, dimensionally, must be in units of feet. 

The Shields critical shear stress equation has been experimentally evaluated in laboratory flumes 
using uniform size particles. Experimental results suggest a constant value in the range of 0.04 
to 0.06 (cited in Reid et al. 1997). Extensive work has been done to evaluate critical shear stress 
in natural channels characterized by a range of conditions (Montgomery and Buffmgton 1997). 

Shear stress assessments typically compare critical shear stress to average total shear stress. 
Average shear stress is calculated by the equation, 

where x is the shear stress, y is the specific weight of water, R is the hydraulic radius, and s is the 
slope of the channel. 

Concerning channel design methods, the engineered channel is typically designed to generate a 
critical shear stress capable of initiating motion in a Dg^ size particle during a bankfull discharge 
event. This concept is consistent with several studies (Pickup 1976; Jackson and Beschta 1992; 
Grant 1987; Carling 1988; Sidle 1998; Booth 1990; and Leopold 1992). Significantly increasing 
or decreasing stream energy by modifying channel dimensions may lead to channel degradation 
or aggradation, respectively. 

C.5.1 Methods 

A critical shear stress analysis was completed for selected reference reaches. Hydraulic models 
were developed using HEC-RAS (US Army Corps of Engineers 2004) to estimate average total 
shear stress. As indicated previously, the model was calibrated by using field-surveyed channel 
geometry, bankfull indicators and bankfull discharges to predict Manning's roughness values. 
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Most available methods for calculating critical shear stress use the Shields equation, but differ in 
their means of calculating critical dimensionless shear stress. There are numerous critical 
dimensionless shear stress equations that have been developed to predict incipient motion. 
Equations are generally derived from laboratory flume experiments that used sand bed channels 
to develop sediment transport relationships. Models based on field data collected in alluvial, 
gravel-bed rivers are less common. The following two methods were selected for their 
applicability to pool-riffle bed forms. 

Parker and Klingeman (1982) 



0.035(D,/D50s) 



■0.94 



Di = Bed particle size for which i% is finer (mm). D84 used. 
Dsos = Bed surface particle size for which 50% is finer (mm) 

Rosgen (unpublished) 

* 

T = 0.0384(Di/D5o)"'^^^^ 

c 

Di = Pavement Dioo or bar sample Dioo (mm) 
D50 = Bed particle size for which 50% is finer 

Allowable boundary shear stress can also be calculated from the following equation developed 
from the Shields diagram (Shields 1936) 



T = pV^/[5.751og(12.27R/ks)]^ 



P 



Density of waster (slugs/ft ) R = Hydraulic radius (ft) 

= Grain roughness (usually taken as 3.5 D84 (ft)) V = Mean velocity (ft/s) 



C.5.2 Results 



Table C-6 presents the results of the analysis. 



Table C-6. Comparison of critical shear stress and actual shear stress (lbs/ft^) for selected reference 
reaches and proposed reaches. 



Shear Stress 
Method 


CFR3 
Middle C4 


CFRl 
Bandmann F3/1 


BFRl 
Bonner B3 


BFR 

Ovando 

C4 


BFR 
Ovando B3 


BFR 
Ovando F4 


Parker and 
Klingeman 

Rosgen 
Shields 


0.63 
0.50 
0.73 


0.82 
1.30 
1.22 


1.57 
0.84 
1.25 


0.77 
0.80 
0.70 


0.97 
0.80 
0.77 


0.67 
1.00 
0.73 


Average of 
Methods (Qbf) 


0.62 


1.12 


1.22 


0.75 


0.85 


0.80 
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Table C-6 (Continued). Comparison of critical shear stress and actual shear stress (lbs/ft^) for selected 
reference reaches and proposed reaches. 


Shear Stress 
Method 


CFR3 
Middle C4 


CFRl 
Bandmann F3/1 


BFRl 
Bonner B3 


BFR 

Ovando 

C4 


BFR 
Ovando B3 


BFR 
Ovando F4 


Actual Shear 
Stress (Qbf) 


0.43 


0.79 


0.95 


0.57 


0.66 


0.62 


Actual Shear 
Stress (Q2) 


0.51 


0.98 


1.16 


0.73 


0.87 


0.79 


Actual Shear 
Stress (Qio) 


0.70 


1.35 


1.57 


0.91 


1.07 


0.96 


Actual Shear 
Stress (Q25) 


0.80 


1.51 


1.73 


1.11 


1.26 


1.13 


Actual Shear 
Stress (Q50) 


0.87 


1.61 


1.84 


1.37 


1.40 


1.32 


Actual Shear 
Stress (Qioo) 


0.92 


1.71 


1.94 


1.65 


1.56 


1.48 


Actual Shear 
Stress (Q500) 


1.04 


1.88 


2.14 


2.46 


1.59 


2.02 


Design Shear 
Stress (Qbf) 


0.53-0.81 


1.19 


0.92 


N/A 


N/A 


N/A 



* Calculated critical shear stress values based on Ds 



C.5.3 Conclusions 



Results yielded an acceptable range of critical shear stress values for bankfull discharge 
conditions. When compared with actual shear stress values, average critical shear stress values 
correspond closely to actual shear stress at Q2 indicating that incipient motion is likely to occur at 
discharges slightly above bankfull. Assuming that the bankfull discharge is responsible for 
incipient motion, it could be concluded that the reference reaches are experiencing aggradation. 
As discussed earlier, significant aggradation was not observed in the reference reaches and 
variation among the results could indicate that dynamic equilibrium can occur within a range of 
shear stress values. In fact, the higher critical shear stress could be a result of bed armoring 
which was observed at all sites. The effects of bed armoring are discussed in more detail in C.7. 
Assuming similar bed composition for the proposed design reaches, design shear stress values at 
Qbf lie within the range of values presented in Table C-4 for average critical shear stress. 

C.6 SEDIMENT TRANSPORT ANALYSIS 

Sediment transport analysis is generally the most challenging component of a hydraulic analysis 
due to the difficulty in collecting bedload data during high flows. Sediment transport data are 
essential for calibrating sediment transport models. Without such data, model calibration may 
not be possible and modeling results should be treated as suspect. 

Despite these difficulties, understanding sediment transport and channel resistance to scour and 
deposition are essential for both understanding existing channel conditions and completing 
channel designs capable of conveying the available flow and sediment load while minimizing 



^ 



River 

DESIGN 

Group . 



C-13 



^ 

i^v 



Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 

channel failure risk. As a design parameter, sediment transport analyses focus on providing for 
sediment continuity, a factor that is necessary for designing a stable channel. In this context, 
channel stability implies that although the channel configuration may change over time, the 
channel's cross-section dimensions remain consistent over time and space. Under the stable 
channel condition, rates of erosion and deposition are approximately balanced as the channel 
erodes its outer banks and stores sediment in bars. 

C.6.1 Methods 

In the absence of bed load data, sediment transport analyses relied on surface and sub surface 
sediment sampling from reference reaches. Sediment data and estimated transport rates from 
upstream and downstream reference reaches formed the boundary conditions for sediment input 
and output rates for the restoration project area. Due to the variability of results from available 
methods, objectives of the analysis focused on demonstrating sediment transport continuity 
rather than predicting actual transport rates. To demonstrate continuity, comparisons were drawn 
between transport rates for upstream reference reaches, proposed project reaches and a 
downstream reference reach. For this analysis, it is assumed that bed material within the project 
reaches will be substantially similar to that found in the respective reference reaches. 

There are numerous equations that have been developed to predict sediment transport rates. 
Equations are generally derived from laboratory flume experiments that used sand bed channels 
to develop sediment transport relationships. Models based on field data collected in alluvial, 
gravel-bed rivers are less common. The following methods were selected for their applicability 
to gravel bed rivers. 

Meyer-Peter and Muller (1948) using D ^ogs 

qb = 8[ps/(Ps-p)] V(g/ p) [(nVnt)''" pSD-0.047(p,-p) Dsoss]''' 

qb = Bed load transport rate per unit width (kg/s) 

Ps = Density of sediment (kg/m'') p = Density of water (kg/m'^) 

g = gravitational acceleration (m/s^) S = Slope (m/m) 

n, = Roughness of bed = S^'^ R^'Vv n' = particle roughness = D9o^'^/26 

R = Hydraulic radius (m) V = Mean velocity (m/s) 

D = Average flow depth (m) 

D90 = Particle size for which 90% is finer (mm) 

D50SS = Subsurface Particle size for which 50% is finer (mm) 

Ackers and White (1973) using Dj 

qb = Ggn f Ps D, V (V/u*)" 8f.[p,/(p,-p)] V(g/ p) [(n'Mt)''' pSD- t * (p.-p) Dsoss]''' 

qb = Bed load transport rate per unit width (kg/s) 

Ps = Density of sediment (kg/m ) 

Ggri = Dimensionless transport rate for i* class 

fi = Proportion of subsurface material in i class (i used = 84) 

Di = Particle size for which i**^ %> is finer (mm) (i used = 84) 

V = Mean velocity (m/s) D = Average flow depth (m) 

u* = Shear velocity = V(gDS) S = Slope (m/m) 



><:^^ DESIGN ( V\r 



C-14 



t^. 



Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam ■ 



October 2005 



C.6.2 Results 

Figures C-3 through C-8 present sediment transport rates for the proposed channels at various 
discharges. Each figure displays a plot of sediment transport rate in relation to valley stationing. 
The intent of this exercise is to evaluate sediment transport continuity between upstream reaches, 
project area reaches and downstream reaches. In general, sediment transport rate is expected to 
increase in the downstream direction as discharge and sediment load increase. 

Station -5000 represents the Bandmann Flats reference reach in the CFR below the confluence. 
Station 5000 represents the design CFR channel below the confluence. Upstream of station 
5000, the BFR and CFR are displayed separately at stations 7500 and 15,000 before entering the 
confluence. At station 7500, the sum of the transport rates for the CFR and BFR approximates 
the required transport rate for the CFR below the confluence at station 5000. Station 18000 
represents an upstream reference reach in CFR3B. 




1,000 



-5000 



5000 10000 15000 20000 

Station (Valley Profile) 



-Ackers & White 



■Meyer-Peter & Muller 



Figure C-3. Sediment Transport Rates for Proposed Channels at Bankfull Discharge. 
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Figure C-4. Sediment Transport Rates for Proposed Channels at Q2 Discharge. 
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Figure C-5. Sediment Transport Rates for Proposed Channels at Qio Discharge. 
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Figure C-6. Sediment Transport Rates for Proposed Channels at Q25 Discharge. 
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Figure C-7. Sediment Transport Rates for Proposed Channels at Q50 Discharge. 
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Figure C-8. Sediment Transport Rates for Proposed Channels at Qioo Discharge. 

C.6.3 Conclusions 

As shown in figures C-3 through C-8 results produced a wide range of values for sediment 
transport rates. As indicated previously, this is common with sediment transport analyses. 
Results were found to be sensitive to slight changes in input variables. Given this sensitivity and 
unacceptable range of values, conclusions focus on degree of sediment continuity between 
upstream and downstream reaches. 

Results indicate that sediment transport capability increases in the downstream direction, 
implying that conveyance is adequate and reaches will not experience aggradation. When 
compared with the upstream reference reach, the proposed CFR3/2 reach appears to be unable to 
transport the potential sediment input from the upstream CFR3 reference at higher discharges. 
Although the difference in transport rates may not be significant, it could be attributed to the 
floodplain constriction that occurs as the CFR3 reach transitions to a more confined valley type 
in CFR2. More detailed modeling and further review of this transition will occur in the final 
design phase. Similarly, the potential sediment output from the proposed CFRl reach appears to 
be greater than the transport capacity of the downstream CFRl Bandmann reference reach at 
higher discharges. This could be due in part to the location of the CFRl Bandmann reference 
below the dam. It is likely that the dam has altered the substrate composition and effects of peak 
discharges on this reach. Again, the difference may be insignificant but it warrants further 
review of potential effects of downstream deposition. 
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Figure C-9. Illustration of sediment transport continuity for selected reference reaches and project 
reaches at bankfull discharge. Values shown are averages values for all methods 

In summary, a range of 3 to 30 percent difference in average sediment transport capability was 
observed. Again, this range is attributed to the sensitivity of the methods to slight changes in 
input variables and transitions in valley morphology. Except for existing conditions, percent 
difference was observed to increase as discharge increased. Results for Qbf and Q2 produced the 
lowest range of differences. Given the sensitivity of the methods and range of differences for 
existing conditions (21 to 23 percent), the analysis appears to provide an adequate demonstration 
of sediment transport continuity between upstream, restoration project area and downstream 
reaches. However, further review of effects of sediment transport discontinuity should be 
evaluated during final design. 
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C.7 SCOUR ANALYSIS 

Several types of scour occur in rivers. General scour is a result of unbalanced sediment 
transport, which may include bed degradation over long periods of time or during extreme flow 
events. Contraction scour occurs when flow accelerates through an abrupt channel or floodplain 
constriction. Local scour occurs in the vicinity of an object that causes turbulent flow such as a 
drop structure, bridge pier, bridge abutment, isolated boulder or fallen tree. Bend scour occurs 
on the outside of channel meanders and is caused by turbulent flow generated by changes in flow 
vectors. 

Scour can occur under live-bed or clear-water conditions. Live-bed conditions exist when the 
sediment transport rate is at or near the sediment transport capacity. Under live-bed conditions, 
erosion is offset by deposition. Clear-water conditions exist when the sediment transport rate is 
significantly less than sediment transport capacity. In such a case, potential exists for erosion 
without deposition. Accordingly, clear-water conditions typically produce greater scour depths 
than live-bed conditions. 

Channel armoring can also affect scour. Armoring is common in gravel bed rivers due to the 
wide range of particle sizes available for transport. Armoring occurs when shear stress is 
insufficient to mobilize a portion of the bed material. Under these conditions, smaller particles 
are entrained and removed from the bed leaving the larger material behind in an armor layer. 
When the armor layer is disturbed as with channel construction, potential exists for scour in the 
reconstructed channel until an armor layer re-forms. 

C.7.1 Methods 

Most scour analysis methods are empirically derived from experiments or field measurements 
and can be biased toward site-specific experiment conditions. Moreover, most methods have 
been calibrated in sand bed channels. Usually, these methods grossly overestimate scour in 
gravel-bed rivers. Efforts focused on using equations that are applicable to gravel-bed rivers. In 
addition, a field-based scour analysis was implemented using scour chains. At the time of this 
report, no scour was recorded due to low peak discharge conditions. 

A scour analysis was completed for the proposed project reaches to evaluate general scour and 
bend scour. General scour was analyzed to determine the required scour depth to re-form the 
armor layer in the new channel. Bend scour was analyzed to validate maximum pool depths. 
The following equation was used to predict the depth of degradation that would need to occur to 
re-form an armor layer. 

U.S. Bureau of Reclamation (1984) 

Ys = Ya(l/Pc-l) 

Ys = Depth of scour required to form armor layer (ft) 

Ya = Thickness of armor layer (ft). Assumed to be twice critical particle size. 

Pc = Fraction of material coarser than the critical particle size 

The following equations were used to predict the depth of bend scour in the design meanders. 
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Thorne(1997) 

Y„,ax = Y„[2.07-0.191og(Rc/Wu-2)] 

Ymax = Maximum water depth in bend (ft) 

Yu = Maximum water depth in crossing upstream of bend (ft) 

Re = Centerhne radius of bend (ft) 

Wu = Water surface width at upstream of bend (ft) 

Mavnord(1996) 

Y„,ax = Y„[l. 8-0.05 l(RcAVu) + 0.0084(Wu/ Yu)] 



U.S. Army Corps of Engineers (1994) 
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Figure C.IO. U.S. Armj' Corps of Engineers momograph for predicting bend scour. 
Table C-7 summarizes the results of the analysis. 
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Table C-7. Summary of predicted scour depths (ft) for selected reference reaches and proposed reaches at 
bankfull discharge. 


Shear Stress 
Method 


CFR3 
Middle C4 


CFRl 
Bandmann 

F3/1 


BFRl 
Bonner B3 


BFR 

Ovando 

C4 


BFR 
Ovando B3 


BFR 
Ovando F4 


Depth of Scour to 
Form Armor 


3.0 


2.8 


1.3 


2.5 


1.1 


2.3 


Depth of Scour to 
Form Armor 
Design Reaches 


2.2 


1.6 


1.3 








Bend Scour - 
Thorne 


5.9-7.9 


11.2- 13.2 


7.7-8.8 


5.9-8.7 


5.3-6.6 


5.2-7.3 


Bend Scour - 
Maynord 


6.3-7.4 


10.3- 11.2 


7.7-8.3 


6.2- 7.2 


5.5-6.0 


5.7-7.3 


Bend Scour - 
USAGE 


7.6-9.9 


11.9-13.5 


8.8 -9.9 


7.9-11.9 


5.2-6.3 


6.1-9.3 


Range of Values 


5.9-9.9 


10.3-13.5 


7.7-9.9 


5.9-11.9 


5.2-6.6 


5.2-9.3 


Avg Bend Scour 


7.5 


11.5 


8.6 


7.4 


5.6 


6.6 


Observed Pool 
Max Depth 


7.5 


23.8 


N/A 


7.8 


7.9 


6.9 


Selected Design 


8.0- 11.1 


18.2-24.3 


9.8- 14.7 


N/A 


N/A 


N/A 



C.7.2 Results 

Results indicate that a significant amount of scour would occur before the armor layer re-forms 
in the design channels. Grade control structures will be required to prevent the reconstructed 
channel bed from significant degradation. Frequency of grade control structures will depend on 
channel slope and threshold for allowable bed degradation. This is discussed in more detail in 
4.5. 

Bend scour results yielded an acceptable range of values for bankfull discharge. Except for the 
CFRl Bandmann reach, average results corresponded closely with observed maximum pool 
depths in the reference reaches. The observed maximum pool depth in the CFRl Bandmann 
reach occurred on the outside of a bend near a bedrock outcrop. A similar bedrock outcrop is 
present in the proposed CFRl design reach near the dam. For this reason, the observed value 
was given greater emphasis when selecting design values. A supplemental check of 100-year 
scour depths was performed for the CFRl and BFRl proposed reaches, yielding values of 24.1 ft 
and 20.3 ft, respectively. Further consideration will be given to 100-yr scour depths and 
associated design implications in Phase III final design. 

C.8 ICE ASSESSMENT 

In Montana, where rivers can develop thick ice covers during the winter, ice jams can contribute 
significantly to channel scour and flooding. Although discharges may be relatively low, the 
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Stages of ice jam flooding may be among the highest on record. Ice jams typically occur 
repeatedly in the same locations and ice jam flooding tends to be localized. 

In the U.S. Army Cold Regions Research and Engineering Laboratory Technical Note entitled 
"Methodology for Ice Jam Analysis" (1980), ice jams were classified as freezeup, breakup, 
floating and grounded jams. Freezeup jams are associated with the formation of frazil ice, which 
eventually forms a continuous ice cover. Breakup jams occur when ice cover is followed by a 
rapid increase in runoff and an ice dam forms. Floating jams are not anchored to the channel bed 
and flow passes beneath the jam. Grounded jams are anchored to the channel bed and flow 
passes over the jam and onto the floodplain. 

Ice jams occur at locations where abrupt changes in ice conveyance occur. Common features 
that contribute to ice jams include transitions from steep to mild slope, confluences, tight 
meanders and channel obstructions. Nearly all of these features are present in the vicinity of the 
restoration project area. The history of ice jams occurrences on the BFR and CFR confirms this. 
Table C.9 presents a summary of historical ice jams in the vicinity of the restoration project area. 



Table C-8. Summary of historical ice jams on BFR and CFR near Milltown Dam.* 




River 


Date 


Location 


Type 


Size 


Ice Stage 


Discharge 
Stage 






y4-mile 






16 ft 




BFR 


02-09-96 


upstream 
from Bonner 


Breakup 


N/A 
8'H 


Max. on 
record 


5.2 ft 


BFR 


02-11-96 


Milltown 


Breakup 


40' W 
4miL 


N/A 


4.3 ft 


CFR 


01-08-30 


Missoula 


N/A 


N/A 


4.41 ft 


2.9 ft 


CFR 


01-11-77 


Missoula 


N/A 


N/A 


Max. 
annual 
7.26 ft 


3.0 ft 


CFR 


11-28-93 


Bonner 


N/A 


N/A 


Max. 
annual 
10.24 ft 


2.4 ft 


CFR 


10-10-95 


Bonner 


N/A 


N/A 


Max. 
annual 
9.05 ft 


3.0 ft 


CFR 


02-08-96 


Turah 


Breakup 


N/A 


Max on 
record 


5.6 ft 



* U.S. Army Research and Development Center Cold Regions Research and Engineering Laboratory Ice 
Jam Database 

Historically, ice jams have formed upstream of Milltown Dam at Stimson Dam and Turah 
Bridge. Upon removal of Stimson Dam, the potential increases for ice jams to move or form 
further downstream. On the BFR, potential areas for ice jam formation include obstructions 
from bridge piers and abutments on the Highway 200 bridge, Burlington Northern Railroad 
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bridge and Interstate 90 bridges. In addition, the confluence could be another potential area for 
ice jam formation. On the CFR, ice jams are likely to reoccur at the Turah bridge. However, ice 
jams that move past the bridge could become lodged in CFR meander bends or steep-to-mild 
transition areas between riffles and pools. The confluence also presents another possible ice jam 
location on the CFR. Records indicate that ice jams on the CFR and BFR account for the 
maximum stage on record. 

EPA has commissioned the U.S. Army Corps of Engineers to perform an ice analysis for pre and 
post dam removal conditions. The results of this analysis are not yet available. During Phase III 
final design, RDG and WWC will rely on the results of this analysis for potential post dam 
removal effects on the stability of the proposed design. 

C.9 HEC-RAS OUTPUT 



River 
Station 



Q Total 
(els) 



MaxClilDpth 
(ft) 



Sliear Clian 
(Ib/sq ft) 



Flow Area Mann Wtd 
(sq ft) ChnI 



TopWChnI 
(ft) 



Hydr Depth C 
(ft) 



VelChnI 
(fl/s) 



Power Chan 
(lb/ft s) 



Bandmann 


660 QBankfull 


10400 


8.57 


0.84 


1455.45 


0.032 


238.21 


6.01 


7.24 


6.06 


Bandmann 


200 QBankfull 


10400 


9.29 


0.81 


1645.71 


0.037 


228.51 


7.15 


6.36 


5.18 


Bandmann 


QBankfull 


10400 


9.21 


0.72 


1311.25 


0.027 


217.18 


6.04 


7.93 


5.68 


Blackfoot Gage 


2228 QBankfull 


6200 


6.49 


1.00 


948.33 


0.038 


194.90 


4.86 


6.44 


6.45 


Blackfoot Gage 


1999 QBankfull 


6200 


7.93 


0.86 


999.27 


0.038 


176.10 


5.86 


6.12 


5.25 


Blackfoot Gage 


1829 QBankfull 


6200 


7.41 


1.28 


962.51 


0.044 


188.70 


5.10 


6.34 


8.11 


Blackfoot Gage 


402 QBankfull 


6200 


10.01 


0.69 


1130.37 


0.038 


195.92 


5.86 


5.48 


3.80 


Blackfoot Gage 


QBankfull 


6200 


6.96 


0.90 


1000.86 


0.038 


200.55 


4.97 


6.12 


5.48 


Ovando B 


550 QBankfull 


3500 


4.77 


0.81 


638.49 


0.038 


162.48 


3.87 


5.56 


4.49 


Ovando B 


QBankfull 


3500 


4.14 


0.51 


725.28 


0.034 


217.31 


3.33 


4.83 


2.47 


Ovando C 


6200 QBankfull 


3500 


6.97 


0.44 


712.58 


0.032 


165.00 


4.26 


4.97 


2.20 


Ovando C 


5340 QBankfull 


3500 


4.76 


0.37 


720.64 


0.028 


187.03 


3.66 


5.09 


1.91 


Ovando C 


250 QBankfull 


3500 


4.93 


0.90 


575.70 


0.036 


167.14 


3.44 


6.08 


5.46 


Ovando C 


QBankfull 


3500 


5.31 


0.56 


595.42 


0.030 


154.31 


3.85 


5.88 


3.32 


Ovando F 


1490 QBankfull 


3500 


5.73 


0.63 


609.55 


0.032 


171.21 


3.56 


5.74 


3.60 


Ovando F 


QBankfull 


3500 


5.33 


0.60 


720.70 


0.038 


174.75 


4.12 


4.86 


2.92 


Ref CFR 3 


1730 QBankfull 


3200 


4.57 


0.28 


696.32 


0.025 


191.94 


3.37 


4.90 


1.39 


Ref CFR 3 


1560 QBankfull 


3200 


4.26 


0.45 


660.76 


0.030 


192.62 


3.24 


5.10 


2.28 


Ref CFR 3 


QBankfull 


3200 


5.14 


0.57 


502.04 


0.027 


142.28 


3.48 


6.46 


3.67 


CFR 3/CFR 2 


14000 QBankfull 


3200 


4.57 


0.5 


545.28 


0.028 


147.88 


3.69 


5.87 


2.91 


CFR 3/CFR 2 


8000 QBankfull 


3200 


4.48 


0.77 


532.28 


0.034 


147.53 


3.61 


6.01 


4.65 


CFR 3/CFR 2 


6000 QBankfull 


3200 


4.62 


0.81 


535.67 


0.035 


146.00 


3.65 


6.00 


4.89 


CFR 1 


1000 QBankfull 


10400 


7.70 


1.18 


1502.56 


0.040 


243.00 


6.18 


6.92 


8.20 


CFR 1 


700 QBankfull 


10400 


7.70 


1.19 


1502.36 


0.040 


243.00 


6.18 


6.92 


8.21 


CFR 1 


400 QBankfull 


10400 


7.69 


1.19 


1501.95 


0.040 


243.00 


6.18 


6.92 


8.21 


BFR 1 


1000 QBankfull 


6200 


6.17 


0.92 


974.06 


0.037 


196.00 


4.97 


6.37 


5.87 


BFR 1 


700 QBankfull 


6200 


6.17 


0.92 


974.27 


0.037 


196.00 


4.97 


6.36 


5.87 


BFR 1 


400 QBankfull 


6200 


6.17 


0.92 


974.14 


0.037 


196.00 


4.97 


6.36 


5.87 
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October 2005 



Reach 


River Profile 


Q Totai 


Max Chi Dpth 


Shear Chan 


Fiow Area 


MannWtd 


Top Wtd 


Hydr Depth C 


Vei Chni 


Power Chan 




Station 


(cfs) 


(ft) 


(ib/sq ft) 


(sq ft) 


Chni 


(tt) 


(tt) 


(ft/s) 


(ib/ft s) 


Bandmann 


660 02 


14900 


10.14 


0.99 


1876.35 


0.032 


277.17 


7.57 


8.19 


8.12 


Bandmann 


200 Q2 


14900 


10.85 


1.05 


2020.97 


0.037 


243.30 


8.71 


7.47 


7.84 


Bandmann 


02 


14900 


10.66 


0.89 


1633.51 


0.027 


225.18 


7.49 


9.15 


8.13 


Blackfoot Gage 


2228 Q2 


8700 


7.77 


1.17 


1217.41 


0.038 


220.32 


6.14 


7.25 


8.52 


Blackfoot Gage 


1999 Q2 


8700 


9.20 


1.09 


1235.81 


0.038 


191.69 


6.93 


7.11 


7.73 


Blackfoot Gage 


1829 02 


8700 


8.71 


1.53 


1216.97 


0.044 


201.36 


5.24 


7.17 


10.97 


Blackfoot Gage 


402 02 


8700 


11.23 


0.88 


1404.76 


0.038 


228.21 


6.85 


6.39 


5.64 


Blackfoot Gaqe 


Q2 


8700 


8.13 


1.11 


1254.19 


0.038 


221.77 


6.13 


7.04 


7.79 


Ovando B 


550 02 


5400 


5.79 


1.08 


888.38 


0.038 


290.11 


4.81 


6.68 


7.23 


Ovando B 


02 


5400 


5.14 


0.66 


944.95 


0.034 


223.69 


4.33 


5.74 


3.80 


Ovando C 


6200 02 


5400 


8.04 


0.62 


906.63 


0.032 


186.59 


5.33 


6.11 


3.80 


Ovando C 


5340 02 


5400 


5.75 


0.49 


1008.88 


0.028 


294.84 


4.62 


6.06 


2,96 


Ovando C 


250 02 


5400 


6.07 


1.09 


776.06 


0.036 


185.84 


4.59 


7.03 


7.68 


Ovando C 


02 


5400 


6.47 


0.72 


789.54 


0.030 


180.21 


4.96 


6.95 


5.03 


Ovando F 


1490 02 


5400 


6.85 


0.79 


808.33 


0.032 


183.01 


4.65 


6.73 


5.30 


Ovando F 


Q2 


5400 


6.56 


0.78 


938.99 


0.038 


180.80 


5.35 


5.78 


4.51 


Ref CFR 3 


1730 02 


4600 


5.42 


0.34 


988.22 


0.025 


402.85 


4.22 


5.53 


1.85 


Ref CFR 3 


1 560 02 


4600 


5.15 


0.50 


930.83 


0.030 


316.01 


4.14 


5.61 


2.80 


Ref CFR 3 


Q2 


4600 


5.97 


0.70 


653.13 


0.027 


201.31 


4.31 


7.44 


5.24 


CFR3/CFR2 


14000 02 


4600 


5.06 


0.55 


1932.04 


0.028 


3001.84 


4.17 


6.29 


3.44 


CFR 3/CFR 2 


8000 02 


4600 


5.25 


0.91 


941.46 


0.034 


612.93 


4.36 


6.74 


6.14 


CFR 3/CFR 2 


6000 02 


4600 


5.43 


0.99 


785.27 


0.035 


316.54 


4.46 


6.86 


6.81 


CFRl 


1 000 Q2 


14900 


9.18 


1.47 


1873.83 


0.040 


255.67 


7.67 


7.99 


11.74 


CFR1 


700 02 


14900 


9.18 


1.47 


1873.80 


0.040 


255.67 


7.67 


7.99 


11.74 


CFRl 


400 02 


14900 


9.18 


1.47 


1873.72 


0.040 


255.66 


7.67 


7.99 


11.74 


BFR1 


1000 02 


8670 


7.28 


1.13 


1194.93 


0.037 


203.06 


6.08 


7.28 


8.21 


BFR1 


700 02 


8670 


7.28 


1.13 


1195.06 


0.037 


203.05 


6.08 


7.28 


8.20 


BFR1 


400 02 


8670 


7.27 


1.13 


1194.73 


0.037 


203.04 


6.07 


7.28 


8.21 



Reach 


River Profile 


Q Totai 


Max Clii Dptii 


Siiear Ciian 


Fiow Area 


i«annWtd 


Top Wtd 


Hydr Deptti C 


Vei Ciini 


Power Ciian 




Station 


(cfs) 


(ft) 


(ib/sq ft) 


(sq ft) 


Ciini 


(ft) 


(ft) 


(ft/s) 


(ib/ft s) 


Bandman 


660 Q10 


25900 


13.31 


1.29 


2840.89 


0.032 


326.81 


10.74 


9.90 


12.77 


Bandman 


200 Q10 


25900 


13.97 


1.54 


2800.56 


0.037 


256.48 


11.83 


9.51 


14.65 


Bandman 


Q10 


25900 


13.58 


1.23 


2354.18 


0.027 


254.91 


10.41 


11.40 


14.05 


Blackfoot Gage 


2228 Q10 


14600 


10.16 


1.50 


1778.33 


0.038 


244.36 


8.54 


8.65 


12.97 


Blackfoot Gage 


1999 Q10 


14600 


11.49 


1.55 


1699.95 


0.038 


214.84 


9.22 


8.92 


13.84 


Blackfoot Gage 


1829 Q10 


14600 


11.04 


2.03 


1721.56 


0.044 


225.14 


8.57 


8.71 


17.68 


Blackfoot Gage 


402 Q10 


14600 


13.50 


1.25 


1931.93 


0.038 


236.57 


9.12 


7.96 


9.93 


Blackfoot Gage 


Q10 


14600 


10.32 


1.50 


1761.29 


0.038 


240.58 


8.33 


8.63 


12.96 


Ovando B 


550 O10 


7800 


6.85 


1.32 


1235.20 


0.038 


351.09 


5.87 


7.63 


10.08 


Ovando B 


Q10 


7800 


6.20 


0.82 


1191.89 


0.034 


242.76 


5.39 


6.65 


5.48 


Ovando C 


6200 O10 


7800 


9.13 


0.84 


1125.11 


0.032 


222.86 


6.42 


7.30 


6.10 


Ovando C 


5340 01 


7800 


6.82 


0.60 


1325.05 


0.028 


294.84 


5.69 


6.95 


4.17 


Ovando C 


250 O10 


7800 


7.23 


1.33 


1006.86 


0.036 


207.28 


5.75 


8.06 


10.73 


Ovando C 


O10 


7800 


7.65 


0.88 


1193.14 


0.030 


563.14 


6.06 


7.94 


7.02 


Ovando F 


1490 O10 


7800 


8.06 


0.95 


1039.40 


0.032 


203.14 


5.87 


7.69 


7.31 


Ovando F 


O10 


7800 


7.87 


0.97 


1179.04 


0.038 


184.48 


6.66 


6.69 


6.51 


Ref CFR 3 


1730 Q10 


9700 


7.88 


0.42 


2113.13 


0.025 


471.12 


6.68 


6.71 


2.84 


Ref CFR 3 


1560 Q10 


9700 


7.69 


0.61 


2012.35 


0.030 


457.90 


6.67 


6.71 


4.09 


Ref CFR 3 


Q10 


9700 


8.23 


1.07 


1211.37 


0.027 


274.53 


6.57 


9.86 


10.58 


CFR 3/CFR 2 


14000 Q10 


9700 


5.97 


0.65 


4658.07 


0.028 


3005.47 


5.08 


7.11 


4.65 


CFR 3/CFR 2 


8000 Q10 


9700 


7.14 


1.12 


2141.19 


0.034 


650.90 


6.26 


7.95 


8.94 


CFR 3/CFR 2 


6000 Q10 


9700 


7.54 


1.46 


1497.53 


0.035 


358.73 


6.57 


8.88 


12.97 


CFRl 


1000 Q10 


25900 


12.16 


2.04 


2671.39 


0.040 


279.51 


10.65 


9.95 


20.31 


CFRl 


700 Q10 


25900 


12.16 


2.04 


2671.07 


0.040 


279.50 


10.65 


9.95 


20.31 


CFRl 


400 Q10 


25900 


12.16 


2.04 


2670.68 


0.040 


279.49 


10.65 


9.95 


20.32 


BFR 1 


1000 Q10 


14600 


9.49 


1.54 


1659.37 


0.037 


216.35 


8.29 


8.95 


13.78 


BFR 1 


700 Q10 


14600 


9.49 


1.54 


1659.50 


0.037 


216.34 


8.29 


8.95 


13.77 


BFR 1 


400 Q10 


14600 


9.49 


1.54 


1659.18 


0.037 


216.33 


8.29 


8.95 


13.78 
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Reach 


River Profile 
Station 


Q Totai 
(cfs) 


Max CM Dpth 
(f1) 


Sliear Clian 
(ib/sq ft) 


Fiow Area 
(sq ft) 


iUlann Wtd 
Clini 


Top Wtd 
(ft) 


Hydr Deptii C 
(ft) 


Vei Chni 

(ft/s) 


Power Clian 
(ib/ft s) 


Bandman 
Bandman 
Bandman 


660 Q25 

200 Q25 

0Q25 


31200 
31200 
31200 


14.63 
15.26 
14.78 


1.41 
1.75 
1.38 


3283.44 
3134.08 
2662.00 


0.032 
0.037 
0.027 


343.38 
261.92 
256.65 


12.06 
13.12 
11.61 


10.54 
10.31 
12.27 


14.83 
18.02 
16.90 



Blackfoot Gage 
Blackfoot Gage 
Blackfoot Gage 
Blackfoot Gage 
Blackfoot Gage 



2228 Q25 
1999 Q25 
1829 Q25 
402 Q25 
Q25 



17300 
17300 
17300 
17300 
17300 



11.12 
12.39 
11.95 
14.40 
11.20 



1.62 
1.74 
2.23 
1.39 
1.66 



2015.45 
1897.33 
1930.16 
2147.16 
1979.42 



0.038 
0.038 
0.044 
0.038 
0.038 



251 .46 
223.87 
230.74 
239.14 
256.66 



9.49 
10.12 
9.49 
10.03 
9.21 



9.16 
9.59 
9.28 
8.55 
9.22 



14.88 
16.68 
20.73 
11.91 
15.30 



Ovando B 
Ovando B 



550 025 
Q25 



10500 
10500 



8.57 



0.038 
0.034 



460.26 
264.54 



6.91 
6.42 



12.96 
7,35 



Ovando C 


6200 025 


10500 


10.21 


1.03 


1366.82 


0.032 


224.50 


Ovando C 


5340 025 


10500 


7.99 


0.68 


1670.41 


0.028 


294.84 


Ovando C 


250 025 


10500 


8.08 


1.72 


1188.03 


0.036 


215.66 


Ovando C 


0O25 


10500 


8.61 


1.01 


1802.56 


0.030 


688.56 



6.59 
6.95 



7.61 
9.39 



5.15 
16.18 



Ovando F 


1490 Q25 


10500 


9.27 


1.10 


1297.24 


0.032 


217.41 


7.07 


8.53 


9.38 


Ovando F 


Q25 


10500 


9.16 


1.16 


1418.22 


0.038 


185.71 


7.95 


7.53 


8.74 


Ref CFR 3 


1730 Q25 


12600 


8.91 


0.48 


2600.22 


0.025 


474.82 


7.71 


7.32 


3.52 


Ref CFR 3 


1560 Q25 


12600 


8.73 


0.68 


2490.26 


0.030 


457.90 


7.72 


7.26 


4.95 


Ref CFR 3 


Q25 


12600 


9.25 


1.24 


1490.13 


0.027 


274.53 


7.58 


10.85 


13.44 


CFR3/CFR2 


14000 Q25 


12600 


6.36 


0.7 


5837.69 


0.028 


3007.04 


5.47 


7.43 


5.17 


CFR3/CFR2 


8000 Q25 


12600 


8.02 


1.2 


2720.48 


0.034 


668.46 


7.14 


840 


10.08 


CFR3/CFR2 


6000 Q25 


12600 


8.48 


1.67 


1845.99 


0.035 


377.66 


7.52 


9.71 


16.25 


CFR1 


1000 Q25 


31200 


13.40 


2.28 


3022.84 


0.040 


289.40 


11.89 


10.70 


24.37 


CFR1 


700 Q25 


31200 


13.40 


2.28 


3022.73 


0.040 


289.39 


11.89 


10.70 


24.37 


CFR1 


400 Q25 


31200 


13.40 


2.28 


3022.55 


0.040 


289.39 


11.89 


10.70 


24.37 


BFR1 


1000 Q25 


17300 


10.37 


1.70 


1851.66 


0.037 


221 .62 


9.17 


9.57 


16.29 


BFR1 


700 Q25 


17300 


10.37 


1.70 


1851.84 


0.037 


221 .61 


9.17 


9.57 


16.28 


BFR1 


400 Q25 


17300 


10.37 


1.70 


1851.59 


0.037 


221 .61 


9.17 


9.57 


16.29 
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Reach 


River Profile 
Station 


Q Totai 
(cfs) 


iUlax cm Dptli 
(ft) 


Siiear Ciian 
(ib/sq ft) 


Fiow Area 
(sqft) 


itiann Wtd 
Ciini 


Top Wtd 
(ft) 


Hydr Deptli C 
(ft) 


Vei Clini 
(ft/s) 


Power Clian 
(ib/ft s) 


Bandman 
Bandman 
Bandman 


660 Q50 

200 Q50 

Q50 


35000 
35000 
35000 


15.53 
16.13 
15.60 


1.48 
1.89 
1.47 


3596.83 
3364.81 
2872.46 


0.032 
0.037 
0.027 


354.64 
270.56 
261.79 


12.96 
13.99 
12.43 


10.95 
10.84 
12.83 


16.22 
20.48 
18.91 


Blackfoot Gage 
Blackfoot Gage 
Blackfoot Gage 
Blackfoot Gage 
Blackfoot Gage 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Appendix C - Preliminary Channel Stability Assessment 

Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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Restoration Plan for the Clark Fork River and Blackfoot River near Milltown Dam - October 2005 
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